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The problem of a spin 1 charged particle with electromagnetic polarizability, 
obeying a generalized 15-component quantum mechanical equation, is investi- 
gated in presence of the external Coulomb potential. With the use of the Wigner's 
functions techniques, separation of variables in the spherical tetrad basis is done 
and the 15-component radial system is given. It is shown that there exists a 
class of quantum states for which the additional characteristics, polarizability, 
does not manifest itself anyhow; at this the energy spectrum of the system coin- 
cides with the known spectrum of the scalar particle. For j = states, a 2-order 
differential equation is derived, it contains an additional potential term r~'^. 

In analogous approach wave functions of the generalized particle are examined 
in presence of external Dirac monopole field. It is shown that there exists one 
special state with minimal conserved quantum number jmin- It this solution, first, 
the polarizability does not exhibits itself, and second, this solution provides us 
with analogue of the known peculiar solution in the theory of the Dirac particle in 
external monopole field, which can describe certain bound quantum mechanical 
state. 

Analysis of the usual vector particle in external Coulomb potential, additional 
to that existing in the literature, is given. It is shown that at j = some bound 
states will arise. The corresponding energy spectrum is found. 



1 Basic equation and notation 



Initial equation for a 15-component wave function of the particle with electric charge and 
polarizability [1-8] has in matrix formalism the form 



(T'^da - m)^ = , ^ 



here blocks of dimensions 1x4, 4x1, 4x6 6x4 respectively are used: 



c 






-pa 































K 










K 














(0) 



(A' 



aN(0) 



5" 



kl 



™ s.k 



{Ka) 

n rib 

The a stands for a free parameter of the model. It concerns with an additional characteristic 
of S" = 1 particle manifesting itself in external electromagnetic and gravitational fields. 

As was shown in [1], equation (1.1) can be extended to general relativity case [9-12] in 
the following way 



[ r"(a;) {da, + B^{x)) - m ] m{x) = 



(2) 
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where 



(3) 



In the present work, some technical possibihties in equations Q underlying by the used 
tetrad formalism will be widely exploited in studying behavior of the generalized vector 
particle in presence of Coulomb and Dirac monopole fields. In that sense, present Section 
helps us prepare for detailed treatment of these two problems. 

Let us consider equation © in Minkowski space assuming the use of spherical coordinates 
and a diagonal tetrad: 



dS^ = c^df-dr^-r\de^ + sm^ 

1 

= (ct, r, e, 0) , ga(3 
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1,0,0,0), ef3) = (0,l,0,0) 
e^i) = (0,0,-,0), ef2) = (1,0,0,^-) 



In tetrad (jU, equation (1.3) will looks 



r°9o + T^dr + 
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where ^ designates 9,4) - dependent operator 
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We will need explicit representation for r"-matrices. With the use of notation 



ei 



1 



62 



1 



63 



1 



1 

















1 


' 63 — 















1 






















1 







-1 





n = 








-1 













, r-i = 


1 













1 







-1 




















the look as 






1 








to 


to 


to 























0* 


0* 





0* 





-/ 





a 




















0* 


0* 





0* 





-/ 





0* 


0* 





0* 


/ 








0* 


0* 





0* 












(4) 



(5) 



(6) 



9. 



n 


n 




n 


to 


to 


to 


n 


n 





n 


to 




to 


u 


u 





u 
















to 





to 


-ej 


to 




0* 





0* 








-n 


0* 


0* 





-el 











0* 


0* 





0* 











Also we will be needing generators of the Lorentz group representations involved: 
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Vector generators are 
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and so on. Tensor generators are (index combinations 01, 02, 03, 23, 31, 12 are used) 
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and so on. Total momentum operators in Cartesian basis has the conventional form 



Jk — h + Sk , Si — iJ^^ , 5*2 — iJ"^^, S; 
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-Sjfc = i [0 © (0 © Tk) © (0 © Tfc) © (Tfe © Tfc) ] 

One should have their representation in spherical tetrad basis. To this end, taking the 
tetrad transformation law 

/„, dx 



where Lorentz transformation is reduced to a pure rotation 



cos 9 cos COS ^ sin — sin 9 

— sin (f) cos (f) 
sin 9 cos sin 9 sin (f) cos 9 



we arrive at the way to connect Cartesian wave function ^ with spherical wave function 

= S ^>{x) , S{9, = 1 © (1 © 0)(1 © O) © (O © O) ] 
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With the help of the same5' one transforms operators Jq = + -S'a, a — 1, 2, 3 : 



Taking in mind the known formulas 

li = i (sin0 de + ctg 9 cos0 9^), 
h — i {— cos 4> de + ctg 9 sin 4> d(t>)-i h — —i d(t> 

and two intermediate ones 

O i de = i T2 , O i d^j, O'^ — i (cos 9 T3 — sin 9 ri] 

for /' = O /„ one derives 
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IL — I2 — i (cos 9 sin 6 T\ — cos T2 + ^ — r- sin 6 T3) , 
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^3 = ^3+ i (sin Ti — cos ^ T3) , 

and additionally for = OraO~^: 

r( = cos 9 cos (p Ti — sin T2 + sin ^ cos , 
T2 = cos 9 sin 4> Ti — cos r2 + sin sin rs , 

= — sin Ti + cos 9 T3 . 

Thus, turns out to be 
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For the following it will be convenient to have diagonal which can be achieved through 
going to the so-called cyclic basis [2]: 
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It is easily verified 
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In this cyclic representation the matrix S3 is diagonal indeed: 

Ss = diag (0; 0, 1, 0, -1; 0, 1, 0, -1; 1, 0, -1; 1, 0, -1 ) . 

Also we will be needing cyclic quantities (all marks reminding about the use of spherical 
tetrad and cyclic representation basis will be omitted below) 

ei = (-^, 0, , 62 = (--^, 0, --^) , ea = (0, 1, 0) , 
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2 Separation of variables for a free particle 

Now we are going to construct proper functions of J^, Ja-operators. In accordance with 
the known differential relations [13] for Wigner's functions 9, 0) 
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the most general form of those proper functions is (see [14]) 
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where the notation is used 

Ds = £>i„,,(0, ^, 0) , s = 0, +1, -1 . 
We will need recurrence relations [13] 

de L>-i = (1/2) ( a L'_2 - ) 
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where v = ^jU + 1) , a = y/{j-l){j + 2) . 

Now let us proceed to separate variables in the main equation (0). Accounting 
structure of all quantities involved after relevant calculation we arrive at 
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Further with the use of intermediate results 
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we get to a radial equation system (the notation v = A/jCT+Ty/v^ is used) 
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Concurrently we will diagonalize P-inversion operator. In Cartesian tetrad it ha s a 
conventional form 



n = [ 1 © (1 © -/) © (1 © -/) © (-/ © +/) ] P , P^(r) = ^(-f) 
which after translating to the spherical-cyclic basis will become 
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The eigenvalue equation n\E' = P \& , with the use of PDg 
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It is easily verified that both (jl7|) and (jl8|) are consistent with the above radial system 
(2.6); at this we will get two sub-systems respectively: 
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3 Vector particle in Coulomb field 

In presence of external electromagnetic fields the original equation Q will become 



[r"(x) {da + Ba{x)-i—Aa) 
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As Coulomb field is described by the potential Aa{x) = (ze/r, 0,0,0) , then most of the 
calculation in previous Sec. 2 and Sec. 3 need not to be repeated - you can immediately 
turn to radial equations with one formal change ( ze^ = a) 
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Concerning eqs. (j^^ for P = (—1)"'+^ not much need to be done to arrive at a final 
second-order differential relation. Actually, from ()22p it follows 
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what is the common equation in the theory of a usual scalar particle in Coulomb field. Its 
solutions are well-known. 

So, in part, the extended theory under consideration behaves like an ordinary 
scalar particle in the Coulomb field. Remaining case when P = (—1)-'+^ turns 
out to be much more involved. For instance, even the simplest case j = leads 
us to an equation too difficult for analytical treatment. 
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Substitution this expression into remaining equation in (j23|) . we get to 
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Unfortunately, we have not been able to proceed with these equations successfully. 

For instance, let us consider in some detail the simplest case of j = 0. For this case a 
special initial form of the wave function must be used (compare it with (jSl ) 
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Operator Tjq^ being applied to that functions vanishes identically; P-operator takes on such 
solutions the proper value —1. So setting z/ = in (|77j) and taking in mind eqs. (pH|) we get 
only three nontrivial equations: 



't + - E2 



, d 2 
- — + - 

dr r 

z(e + -)E2 + ta-^ (^^ 
r dr m 



+ -) E2 = m^o , 



= 



E2) = m $2 , = 



«N I d ^ „ 
- e + - $2 - — $0 = m ^2 , 
r dr 

Excluding variables $0? '^'2, we arrive at an equation for E2: 



= 



(29) 



E2 = - e + - 
r 



..a, .da 

t{e + -) + i<y ^ 

r dr 



Eo 



a a , a. ^ , d 2, ^ 



from which it follows a 2-oder differential equation 
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4 Particle in magnetic monopole field 

Now we will consider the vector particle in external Dirac monopole field [15]. Initial 
wave equation is 
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To describe monopole field the Schwinger monopole potential will be used; the latter after 
translating to spherical coordinates will have the explicit form = (0,0,0, = gcosO) ; 
g stands for a magnetic charge. According to this you are to make one substitution 
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You may immediately pass to equations of the type (fTUI) : 
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In ? 1 / ds + (t'} + ik) cos 9 . ^ -> , , 

r r sm fc' 

It is the point to establish a suitable form of the wave function being proper one of 
rotation symmetry-based operators. In presence monopole background there exist three 
relevant ones: 
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So, the most general form of 15-component wave function with quantum numbers e, j, m will 
look as 
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Now you should substitute the wave function into (IHHj) and exclude variable t,9, 
this you will need some known recurrence relations [13] 
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- n 9e + T2i-t) . n H = m^ 

r sine 

1 r ^+ „ ^+ . .X — — k cos e 



- [ el de + 



^n = m E 



eg* dr $0 

r sm e 

In ? 1 r / -x — — (A; — s-i) cos e -, -? - 
r r sm u 



where 



S3 



+ 10 

0-1 



(36) 



At 



(37) 



(3J 
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with the use of several intermediate formulas 



deC = e 



-iet 



CideDk-i 




C2deDk 









-m- (A; -S3) cos g ^ ^ ^~iet 



sini 



Ci{a Dk-2 - c Dk) 
C2{c Dk-i -dDk+i) 

C,{d Dk-bDk+2) 
Ci{-a Dk-2 -cDk) 
C2(-c Dk-i - d Dk+i] 
Cs{-d Dk-hDk+2) 



e^deC = e-'^' -j= [ -Ci(a Dk-2 - c Du) + {dD^-h Dk+2) ] 

— m — (A; — S3) cos^^ ^ 
-^e2 6 



sin 9 



-iet 



= -j^[-Ci {-a Dk-2 -cDk) - C3 i-d Dk - b Dk+2) ] 
(el9.-.e,^^^^^^^:i^^)C = e--*y2(cC, + dC,)Dk 



sin 9 

r^ / N — m — (A; — S3) cos6' -L ■., . /-. 

ndeH + r2 i-i) . n ^ = e""* ^^2 

sm 9 



(9e + e2*(-z) — )C = e V2 



sin 9 



from ()38|) you arrive at 



-c i/2 /^fc-l 
cH^-dH,) Dk 
+d H2 Dk+i 

-c C(r) Dk-i 


-rf C(r) Dk+i 



- le Co - + -) C2 - — {c Ci + d C3) = m C 
dr r r 



-{^ + -)E2- — icE, + dE,) = mCo 
dr r r 

+«e El + 1{— + -) Hi + H2 = mCi 

dr r r 



+ieE2 



iV2 



[c Hi — d H^) = m C2 



^ , d 1. 
+te E-i - z — + - H^ 
dr r 



iV2d 



Ho = m C3 



-iea C - + -) E2- — (c El + d E3) = m <!>o ■ 
dr r r 

+2e El + 1{— + -) Hi-a C + H2 = m^i, 

dr r r r 

+ie E2 + (T - ^^(c Hi - d Hi) = m ^2 ■ 

dr r 

■ r. •/ ^ ^ rr V2 d ^ iV2 d 

+ie Es - z — + -) Hi - a C H2 = m 

dr r r r 



v^c d 
ie $1 H ^o^m El , -ie $2 — r ^2 

r ar 



E2 , 



-ze $3 H $0 = , -«(^ + -) $1 ^2^rnHi , 

r dr r r 

c $1 - d $3 = m , +«^r + - *3 + ^2^mE^. 42 

r or r r 

5 States with j = j/Vnin in monopole field 

As known, in the monopole problem there arises a very peculiar situation for states with 
minimal value of j. Now we are going to consider those in our case. In accordance with 
quantization conditions by Dirac-Schwinger the parameter k — eg /he involved will take on 
the values 

A; = ±l,±3/2,±2,... and j =| A; | -1, | A; |, | A; | +1, ... (43) 
One should consider with more details several particular cases. Firstly, let k be 

^ — jmin — ■ 

Correspondingly, a wave function should be taken in the form 

: C{t) = , Co(r) = , $o(r) = , 



C{x) = e 
E{x) = e 



—let 



-itt 



Ci{r) 



Ei{r) 





—iet 



, $(x) = e 
, H{x) = e-''' 



$i(r) 



Hi{r) 





(44) 



On those ^^~^{x) the angular operator vanishes identically E^_^^'-'~°(x) = , and corre- 
sponding radial equations are as follows 

+ ie Ei + i{-^ + -) Hi = mCi , 
dr r 

+ieEi+i{-^ + -) Hi^m^i , 
dr r 

d 1 

-i e $1 = m £^1 , -i (— + -) $1 = m i?i . 

dr r 



Excluding i/i, £"1 for Ci = $1 one gets 



. d'^ Ox $1 
dr" 







$1 



exp [ iVm^ — 



(45) 



One of these solutions might be associated with a "bound state", this is what makes them 
very interesting from physical standpoint. 
Another variant 



k — 1, jmin — 
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will look the same. Now an initial substitution is 

.1=0 



: C(r) = 0, Co(r) = 0, <l>o(r)=0 




Csir) 






C{x) = e 



E{x) 



-let 



-let 



H{x) 



-iet 



-iet 













(46) 



Again on those functions {x), ^ (x) = , and radial system will be 

+ i e ^3 - + H3 = mCs, 

ar r 

ar r 

-i e (l>3 = m Es , i (— + -) $3 = m ifg . 

ar r 



which after excluding H^, E gives 



— = exp[ ±vm^ — 



(47) 



Now we consider states with minimal jmin =| ^ I ^1 when k = ±3/2, ±2, ±5/2, ... 
Here, though the 9, 0- variables enters v]/-?™™ (a;) , however the angular operator acts on these 
functions as zero-like. Actually, as 



A; = ±3/2, ±2, ±5/2,... , 



k 



^^™'»(x) : C{r) = 0, Co{r) = 0, %{r) = 



iet 



C{x) = e 
E{x) = e-''' 



Ci{r)Dk^, 



Ei{r)Dk-i 





$i(r)A-i 



Hi{r)Dk-i 





(4J 



Here you need the recurrence formulas [2] 



k~2 



m — [k — 1) cos 6 



sin 9 



Dk-i 



k-1 



-D 



k-2 



with the help of which you easy can prove T^g^^ vl/i™i"(a;) = and then again we arrive at 
eqs. (jl^ . 

Another variant 



k = -3/2,-2,-5/2 



! ••• ! Jmm 



k I -1 



looks the same: 



(x) : C(r) = 0, Co(r) = 0, <l>o(r) = 



Cix) = e 



-let 



E(x) = e 



-iet 












= e' 



H(x) = e 



-iet 



-iet 









i^3(r)A+i 



(49) 



— m — (A; + 1) COS0 



D 



k+2 



sin^ 



Dk+i 



ik + 1] 



D 



A:+2 , 



with the help of which you easy can prove S^,,/, (x) = and then again we arrive at 
eqs. (jl7|). 



6 On the Coulomb problem of an ordinary S = 1 par- 
ticle 

Quantum mechanical study of a vector particle in Coulomb field has a long history. We are 
going to review this old problem else one time in the frame of the approach developed above. 
So, let us consider a 10-component Duffin-Kemmer equation for a S=l particle in Coulomn 
field all details on the working technic used here and based on the tetrad formalism). 

In the spherical-cyclic basis the most general form of 10-component wave function with 
quantum numbers e, j, m is 

*(x) = {$o(x), $(x), E{x),H{x)] , 

$i(r) /^-i 



<l>o(x) = e-^^* $0(0;) Do , $(x) 



E 



-iet 



E2ir) Do 



H 



-iet 



— iet 



$2(r) Do 
$3(r) Z^+i 

Hi{r) 
H2{t) Do 



(50) 



where D^ = Di 



0,0), s = 0, +1, —1 . Omitting details of separation of variables let 



us immediately write down a radial system ( that can be derived directly from 15-component 
above by setting a = and so on) 

- (-^ + -) ^2 - - (^1 + ^3) = m% , 
dr r r 

«N ^ , d 1. 

+t{e + -)Ei +1 {— + -) Hi + i- H2=m^i , 
r dr r r 

a, V , 

+i{e+-)E2 - i- {Hi - H3) = m<^2 , 

+z(e + -)E3 -ti-^ + -)H;-i-H2 = m<l?3, (51) 
r dr r r 



Ifi 



(X u 

-i(e + -)$i +-$0-^^1=0, 

Q: d 

-lie + -) $2 r^o- mE2 = , 

r dr 

-lie + -) $3 + - $0 - rnE^ = , 
r r 

-i (4- + -)*!- i- $2 - mi/i = , 
dr r r 

+i- ($1 - $3) - mi/2 = , 

r 

+i (4- + -) $3 + - mH^ = . 
dr r r 



(52) 



Concurrently with J^^Jz one may diagonalize the operator of the space inversion H. That 
in chosen representation has the form 



n = [ (1 e Hg) e {n^ e -n3) ] p , n3 



0-1 

0-1 
-100 



The eigenvalue equation fl^ — P ^ gives solutions of two types: 

P = {-iy^\ $0 = 0, $3 = -$i , $2 = 0, 
E3 = —El , £'2 = 0, Hs = Hi ; 

P = (-1)^' , E3 = +Ei , Hs^-Hi, H2^0. 

As for the case P — (— l)^"*"-*^ you have the system 

«N ^ , d 1, u 
+ z e + -) Ei+ t{— + -)Hi + 1-H2 = m$i . 
r dr r r 

-lie +-)$! = mEi , -i{4- + -)$i = "^ifi , 2i-$i = mif2 

r dr r r 



(53) 



which after excluding Ei,Hi, H2 leads to 

d^ 2d 



dr^ r (ir r 



$1=0 



(54) 



This exactly the same equation that arises in the theory of a scalar particle in Coulomb field. 

Its energy spectrum and wave functions are well known. 
States with P = {—ly are characterized by the system 

(4- + -)E2 - 2-Ei = m$o , 
dr r r 

+i{e + -) Ei + + -)Hi = m$i , 
r dr r 

a, V 
+i(e + -)E2 - 2%-Ei = m$2 , 

-lie + -)$! + -$o = fnEi , 
r r 

, OLs. ^ d ^ ^ 
-lie + -)$2 - ^$0 = rnE2 , 
r dr 

-i{4- + - ^ mHi (55) 
dr r r 
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in solving which we have not been able to succeed. 

The above radial equations as well as substitutions (j^UI) are correct only for j = 
1,2,... However you should consider the j = case separately and having started with a 
special wave function: 



E(x) 



'let 



-let 





Co(r) 



E2{r) 




H[x) = e-''' 





$2(r) 




H2{r) 




(56) 



This \1/ is eigenfunction of fl with 11 = +1. Corresponding radial system is ($o = V'o, —i^i 



H. = 



a 



e H — )E2 = m(p2 



d 2 

-(— + -)E2 = rriLpo , 
ar r 

a d 
[e + -)lp2 - — v?o = mE2 
r dr 



from where it follows (that corresponds to (|3(J|) when (7 = 0) 



rf2 2d 



+ ^+ e+- 



dr'^ r dr 
or introducing E2{r) = r~'^f{r) 



m 



Eo = 



/ + ( - + — ) / = 



dr 

The latter will read in usual dimension units as 



(57) 



dr 



2f + 



E^ M^c^ ^ 2Z{eych){E/ch) _ 2-Z^{e^/diy 



f = 



c^h^ r r^ 

Now you can pass to a dimensionless coordinate x = rE/ch, then the equation becomes 



d^ 



dx"^' 

With the notation 



1 - ^ + 2Z{eych) 



1 2 - Z^eychf 



x 



x^ 



fix) = 



1 

ch 137 ' 



the equation will read as 



dx^ X 
To obtain bound states you should take f{x) in the form 

fix) = e"^" F{x) 



(5^ 
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where a and h are expected to be positive whereas F{x) be a polynomial in x. From ()58|) we 
derive 

X F" + (2a - 2hx) F' + 
+ [ ^ '—^ + (6^ + 1 - A^)x + (27 - 2a6) ]F = 0. 

X 

This equation with demands 

a(a - 1) + 7^ - 2 = , 6^ + 1 - = 

leads to 

X F" + 2(a - bx) F' + 2(7 - a6) F = . (5.15c) 

For a and 6 we have 



a = l±^E^, . = ±v/A^=±^l!^. (59) 
2 E 

The choice of "upper" signs is what you need to a and b be positive. 
Further, taking F{x) as a series 

00 00 00 

F(a;) = ^ Cfc x\ i^' = A;Cfc x^"^ , = ^ A;(A; - x'^'^ , 

fc=0 k=l k=2 

you get 

00 00 00 00 

J2 Kk - x''-^ + 2a kCk x''~^ - 26 '^^fc 2;'^ + 2(7 - ab) ^ x'^ = 

k=2 k=l k=l k=0 

or 

n{n + l)C„+i x" + 2a ^ (n + l)C„+i x" - 

n=l n=0 

00 00 

- 26 J] nCn x" + 2(7 - ab) ^ C„ x" = 



71=1 n=0 



or 



[ 2aCi + 2(7 - ab) ]x° + [2C2 + 2a 2C2 -2bCi + 2(7 - ab)Ci ] x + 

00 

+ Y[n{n + l)C„+i + 2a(n + l)C„+i - 26nC„ + 2(7 - a6)C„ ] x" = . 



n=2 



Demanding coefficients at all x be equal to zero, we produce recursive relations 

Ci = -(7-a6)Co = 0, 
C2 2 (1 + 2a) = 2 [6 - (7 - ah)] Ci = 0, n = 2, 3, 4, 
C^+i (n + 1) (n + 2a) = 2 [n 6 - (7 - a6) ] = . (60) 

To terminate an infinite series at a certain place you need impose 

Cn+1 = ^ [ 6 - (7 - a6) ] = 



1Q 



that provides us with a quantization condition to produce an energy spectrum. It looks as 



7 — 06 



= N 



which having remembered a and b (5.16) will take the form (the notation 2T = (1 + ^/Q — 47^ 
is used) , 



276 — Ty/m'^c^ — 



2Vr 



,2M _ ^2 



from which you arrive at the energy formula 



e = mc 



1 + 



r 



(r + Ny 



N 



-1/2 



(61) 



Some additional remarks might be given. The material of this Section might be of some 
interest in the light of the well-known old results obtained by I.E. Tamm [...] on behavior of 
the vector particle wave functions in presence of external Coulomb field. 

The most principal part of this work consists in the following. All possible solutions can 
be divided into two classes. One is reduced to the well known functions of the scalar particle 
in Coulomb field with usual and expected behavior at the origin r = and at infinity r = cx). 
However another part of wave functions provides us with big surprise - they arc singular 
at the origin, corresponding differential equations have not been resolved up to now; it is 
generally accepted that they should be associated with a situation of falling down to centre 
r = 0. 

In essence, results of this Section prove that there exists definite correlation between 
values of j, of P, and property of wave function to be singular or non-singular; it can be 
seen in the following tabic: 



Value 


p = i-iy 


p = i-iy 


P = (-l)J+i 


P = (-l)J+i 


of quantum number 


singular 


non-singular 


cingular 


non-singular 




solution 


solution 


solution 


solution 


j = 




P = +l 






i = i 


P = -1 






P = +l 


J = 2 


P = +l 






P = -1 


J = 3 


P = -1 






P = +1 


J =4 


P = +l 






P = -1 
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